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ABSTRACT 

What is the mass of the progenitor of the Sagittarius (Sgr) dwarf galaxy? Here, we reassemble the 
stellar debris using SDSS and 2MASS data to find the total luminosity and likely mass. We find that 
the luminosity is in the range 9.6 — 13.2 x 1Q 7 Lq or My ~ —15.1 — 15.5, with 70% of the light residing 
in the debris streams. The progenitor is somewhat fainter than the present-day Small Magellanic 
Cloud, and comparable in brightness to the M31 dwarf spheroidals NGC 147 and NGC 185. Using 
cosmologically motivated models, we estimate that the mass of Sgr's dark matter halo prior to tidal 
disruption was ~ 1O 1O M0. 

Subject headings: dwarf galaxies: individual (Sagittarius) - dwarf galaxies: luminosity - dwarf galaxies: 
mass 
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1. INTRODUCTION 

The Sagittarius (Sgr) dwarf galaxy is the closest 
Milky Way companion (~ 24 kpc from the sun) and 
is currently being disrupted in the Galaxy's tidal field. 
Both leading and trailing tidal streams have been de- 
tected using many different tr acer populatio ns, includ- 
ing A stars (lYannv et al.ll2000ft , RR Lyraes (I Vivas et al 



2001 1 : iWatkins et al.ll2009ft. r ed clump stars (jMateo et al 
1998; Maicws ki et al.1 11999ft and main sequence turn- 



off st ars ( Martmez-Delgado et al.l 120011: iNewberg et al.l 

2002) . Perhaps the most spectacular views of the Sgr 
stream are given by the M giants (Majewski et al. 

2003) . A portion of the stream around the North 
Galactic Cap has been traced to much fainter mag- 
nitudes using data from the Sloan Digital Sky Sur- 
vey (SDSS, see Belokurov et al. 2006a "The Field of 
Streams") revealing a prominent bifurcation into the 
so-called streams A and B. There has been consider- 
able effort in modelling the disruption of Sgr since the 
stream stars are useful tracers for the shape of the Milky 
Way's dark halo potential. Nonetheless, this work has 
yielded ambiguo us results. Depending on which d atasets 
are fi t, oblate (I Johnston et al.l 2005ft. prolate ([Helmil 



12004ft. spherical (|Fellhauer et al 



2006ft . or even triax- 



ial (jLaw et al.l[2009h halos are preferred. 

One reason for this ambiguity is that the structure, 
size and nature of the Sgr progenitor are very uncertain. 
The original luminosity and mass of Sgr is larger than its 
luminosity and mass today by amounts depending on the 
mass density (including any dark matter) and orbital his- 
tory of the system. A neat illustration of this is provided 
by the work of Jiang & Binney (2000), who used N-body 
simulations to show that the current configuration of the 
core of Sgr can be produced by a very wide variety of 
initial conditions. At one extreme, the Sgr dwarf might 
initially possess a mass in excess of ~ 10 and fall 

into the Galaxy from a distance in excess of 200 kpc. At 
the other extreme, it might be only ~ 10 9 M Q and start 
off at distances similar to its present apocenter of 60 kpc. 

In this paper, we analyse the number of stars and the 
luminosity of the present-day Sgr core and debris. By re- 



assembling Sgr, we provide new constraints on its param- 
eters prior to disruption. We describe the data sets used 
for this study in §2 and analyse the colour-magnitude 
distributions of the stream and core as well as the lumi- 
nosity density along the leading and trailing arms in §3. 
We discuss the implications of our work for estimating 
the mass of Sagittarius and conclude in §4. 

2. DATA 

We trace the Sgr stream using main sequence, red giant 
and horizontal branch stars from the SDSS (York et al. 
2000), together with M giants in the Two Micron All-Sky 
Survey (2MASS). SDSS is an imaging and spectroscopic 
survey that covers one-quarter of the celestial sphere in- 
cluding an area of ~ 8400 deg 2 around the North Galac- 
tic Cap, which contains a portion of the leading arm of 
the Sgr stream. There are also some imaging stripes 
that sparsely sample low Galactic latitude fields (~ 3200 
deg 2 ) which contain parts of the trailing arm. The imag- 
ing data are collected in five passba nds (u,q,r,i,z), for 
which there are model isochrones (Girardi et al.1 120041 : 
iDotter et al.l 12008ft and photometric transformations Q 
readily available. Here, we use positions and magnitudes 
of sta rs extracted from Data Release 7 (jAbazajian et al.l 
2009ft and corrected for extinction using the maps of 
Schlegel. Finkbeiner. fc Davis! ([1998). The bright limit 



in the it-band is about 14 mag, whilst the faint limit is 
about 22. The final SDSS sample contains ~ 3,740,000 
stars along the Sagittarius stream and ~ 3, 200, 000 stars 
in a comparison field in the North Galactic Cap region as 
well as ~ 294, 000 stars along the stream and ~ 255, 000 
stars in the comparison fields in the southern stripes. We 
additionally use the positions and J, H , K magnitudes of 
~ 1,450,000 stars in the 2MASS All-Sky Point Source 
Catalog. This is complete down t o roughly K = 15 mag . 

Finally, we use a dataset from Bcllazzi ni et al.l (2006) 
to study the properties of the Sgr core. It comprises of 
two fields (~ 299, 000 stars on and ~ 57, 000 stars off the 
Sgr core) in two passbands B and V. The sample goes 
approximately 3 magnitudes below the main sequence 
turn-off and includes stars up to the red clump, blue 



mno , vasily, nwe , j or pcgaOast . cam .ac.uk 



1 http:/ /www. sdss.org/dr4/algorithms/sdssUBVRITransform. html 
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horizontal branch and M giants. 

3. ALL-SKY VIEWS OF THE SGR STREAM 

We begin by giving an overview of different observa- 
tions of Sgr and the terminology that we will employ 
throughout the paper. We refer to the still bound por- 
tion of the dwarf as the core or remnant, which we define 
as stars within 10 kpc of the center. These stars follow 
a King profile (see e.g., Majewski et al. 2003). The lu- 
minosity of Sgr is usually quoted as the light integrated 
inside this radius, but it is not obvious that all of these 
stars are still gravitationally bound. The tidal debris 
from the disruption wraps around the Galaxy (possibly 
multiple times) in what is usually referred to as the lead- 
ing and trailing tidal streams or arms. The portion of 
the debris that is seen around the North Galactic Cap is 
bifurcated into stream A and stream B. The working hy- 
pothesis for this paper is that both stream A and stream 
B are part of the leading arm. Part of the leading arm is 
also detected in the 2MASS survey. The trailing arm is 
in part covered by the southern stripes of the SDSS and 
more fully covered in 2MASS. 

The best panorama of the Sgr stream is given by 
M giants selec t ed fro m 2MASS, as first realized by 
IMaiewski et al.l (|2003[) . Their choice of cuts on J, H, K 
magnitudes picks out bright red stars at the typical dis- 
tances of Sgr. For Figure[TJ and in the following analysis, 
we adopt the selection criteria introduced by Majewski 
et al. (2003), namely 

J-K > 0.85 
J -H < 0.561(J-iT) + 0.36 
J-H > 0.561(J-#) + 0.22 (1) 
|6| > 10° 
E(B - V) < 0.555 

Additionally, for the left panel of Figure Q] only, we re- 
move some of the reddest stars with J — K > 1.25 
and limit the magnitude (heliocentric distance) range to 
10 < K < 12 (7 < D < 70 kpc) for the trailing arm 
and to 11 < K < 13 (5 < D < 110 kpc) for the lead- 
ing arm. This is in order to isolate the tidal debris of 
the Sgr. Notice too that the Magellanic Clouds are very 
prominent. Superficially, Sgr already looks to be inter- 
mediate between the Large and Small Cloud in size and 
luminosity. 

The right panel of Figure [1] shows a zoom-in of the 
Sgr core, made using the cuts of eqn (TTJ), together with 
an additional color-magnitude (CMD) box. This is de- 
rived by shifting the ridgeline of Sgr M giants (K — 
-8.650(J - K) + 20.374, Equation 5 of Majewski et 
al. 2003) by ±3.5er or ±1.26 mag in the color range 
0.95 < J — K < 1.1. Using this selection gives a clear 
view of the Sgr core, the photometric center of which is 
m arked. Also shown are the on-core and off-core fields 
of iBellazzini et al.l (|2006[) . 

The SDSS view of the Sgr stream is confined primarily 
to the North Galactic Cap and hence is shown in equato- 
rial coordinates in Figure [2] to minimize dist ortion. The 
"Field of Streams" (jBelokurov et al.l l2006a|) shows the 
Sgr stream over a 100° arc, with an unprecedented level 
of detail. The datasets provided by 2MASS and SDSS 
are complementary - the 2MASS survey picks up the 



trailing arms with good clarity, but it loses the leading 
arm around the North Galactic Cap. 

3.1. Distances 

There is information about the dista nces to debris at 
a number of points along the stream (iBelokurov et al.1 
l2006atlNewberg et al.ll2007l: IWatkins et al.ll2009fl . Helio- 
ccntric di stances to the trailing arm were originally com- 
puted by IMaiewski et all (|2003[ ). who showed that there 
was little distance gradient along this part of the stream. 
There has been more controversy about distances to the 
leading arm, mainly because there are much steeper dis- 
tance gradients. 

Given the crucial role that distances play, we begin by 
seeking independent confirmation of earlier results. For 
the leading arm in the Field of Streams, we select blue 
horizontal branch (BHB) candidates from the SDSS in 
the on and off-stream fie l ds us ing the criteria given in 
Figure 10 of lYannv et all pOOOl) . This gives us ~ 18700 
candidates in the Sgr stream and ~ 8500 candidates in 
the comparison field. We note that these numbers do not 
include all genuine BHBs, and suffer from some contam- 
ination from blue stragglers (BS). The upper panels of 
Figure [3] show the density of BHBs in the two fields, to- 
gether with their difference. The distances to BHBs are 
calculated assuming an absolute magnitude of M g = 0.7, 
which is appropriate for our color range (u — g > 0.9, see 
Table 2 o flSirko et al.l (|2004h ). Running diagonally across 
the plot is the signal from the BHBs in the Sgr stream. 
It is accompanied by a fuzzy shadow of BSs two magni- 
tudes fainter. The Galactic halo foreground appears to 
be slightly asymmetric, most likely due to the presence of 
the Virgo Overdensity (jJuric et al.ll2008| ) in parts of the 
stream field. However, the subtraction works well and 
leaves the Sgr stream clearly visible. The stream signal is 
roughly linear in the plane of distance versus right ascen- 
sion and can be described by a straight line with gradient 
0.386. The solid lines in the two leftmost panels show 
the straight line fit offset by ±5 kpc. The lower panel 
compares the BHB distances (straight l ine) to the lead- 
ing arm with those originally inferred bvfBclokuro v et al.1 
( 2006a) using the i magnitude of the subgiant branch in 
the "Field of Streams" . It is reassuring that there is 
excellent agreement between the two independently cal- 
culated distances, especially as the BHB's luminosity is 
almost independent of age and metallicity. 

Also shown are the distances to 2MASS M giants in 
both leading and trailing arms. Using the cuts described 
above, we create a density distribution of M giant stars 
in Sgr's orbital plane as shown in Figure [4] Distances 
to pixels with highest density along the stream are then 
computed with a typical uncertainty of 6 kpc. This is 
a combination of the uncertainty in the color-absolute 
magnitude relation for M giants and the intrinsic width 
of the stream. Finally, we also show two detections of 
the trailing debris in the SDSS southern stripes. The 
distances are computed using both BHBs and main- 
sequence turn-off stars. The overlaps of the SDSS and 
2MASS datasets will play an important role in what fol- 
lows. 

As a useful summary of the data, we list distances and 
right ascensions in Table[Q Recently. lYannv et alj (|2009Q 
reviewed the detections of the Sgr stream in SDSS data, 
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Fig. 1. — Left: M giants selected from 2MASS tracing the Sgr stream in Galactic coordinates. The box marks the central parts of the 
Sgr dwarf and the arrow indicates the direction of motion of the dwarf. The leading and trailing arms are visible at positive and negative 
Galactic latitud e respe ctively. The gray-scale shows the SDSS footprint. Right: Zoom-in of the area outlined by the box, with the fields of 
Bcllazzini ct al. (2006) marked with squares. The photometric center of Sgr is marked with a star. 
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Fig. 2. — "Field of Streams" (Bclokurov ct al. 2006a) with the Sgr debris bifurcated into a lower (A) and an higher (B) declination 
stream. Selection boxes along the stream (used in §3.3) are shown, together with comparison boxes at the same Galactic latitude but 
reversing the sign in longitude, (i.e., effectively reflecting the boxes in I = 180°). The stream boxes are 5° wide in right ascension and are 
chosen to cover the entire width of the stream in declination. 
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Fig. 3. — Upper: Density of SDSS BHBs in the plane of distance versus right ascension for the Sgr stream (left), comparison field (middle) 
and their difference (right). Note the prominent linear feature associated with streams A and B, which is shadowed by the blue stragglers 
some 2 magnitudes fainter. The solid lines bracket the distance gradient in Sgr's BHBs. Lower: Comparison of distances to the Sgr stream 
in 2MASS (open symbols) and SDSS (filled symbols, from Belokurov ct al. 2006) datasets with the uncertainty in distance shown in 
gray. The leading arm is marked with circles, the trailing with triangles. Note the good match between distances derived from different 
indicators. Arrows mark locations of data-set overlaps. 
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Fig. 4. — Density distribution of 2MASS M giants in the orbital 
plane of Sgr as defined in Majewski et al. (2003, see Figure 11). 
The center of Sgr is marked with a star. The solid red line marks 
the ridge line, whilst the dashed lines indicate the 3<r boundary 
of the tidal arms. The circles are the SDSS observations around 
the North Galactic Cap. The triangles mark the stream detections 
in the SDSS southern stripes. The inset shows the number of M 
giants as a function of heliocentric distance in the dot-dashed slice 
of the stream. A Gaussian fit to the distribution yields a width of 
6 kpc. The Galactic Center and the Sun positions are marked. 

and their results are in good agreement with ours. 

3.2. Sagittarius Core 

In order to derive the luminosity of Sgr, we need to dis- 
entangle its stars from the Galactic fore and background. 
We start by re-visiting the calculatio n of the lum i nosity 
of the remnant. The original value of llbata et al.l (1995) 
is based on data from Schmidt photographic plates. With 
an estimate of the extent of the remnant, together with 
a CMD of a 2.5° x 2.5° central field, they calculated 
an absolute magnitude of My ~ —13. This relies on 
the similarity of the CMDs of Sgr o n the one ha n d and 
Fornax and the SMC on the other. IMateo et all (|1998n 
acquired data on 24 fields of size 0.25° x 0.25° across 
the face of Sgr and built a number density profile of 
main sequence stars . Then, using a su rface brightness 
normalization from IMateo et al.l (|1995f) . this can be in- 
tegrated to give a total luminosity of the remnant as 
My — —13.3. However, due to the s i ze of the field of 
view (0.25° x 0.25°) in IMateo et al.l ([l995l ). the lumi- 
nosity function (LF) does not extend brighter than the 
red clump and does not include BHBs and BSs. So, the 
shape of the surf ace brightness pr o file is well-constrained, 
as confirmed by iMaiewski et al.l (|2003l ) , but the overall 
normalization is more uncerta in. 

The data of IBellazzini et all (12 006) cover a larger area 
than those of IMateo et al.l (|1995l) at the same location, 
and hence give a clearer view of all the stellar pop- 
ulations. The core data are obtained in a 1° x 1° 
field approximately 2° east of the Sgr center located at 



(I, b) ?s (6.5°, —16.5°) and along the major axis. The con- 
trol field is 0.5° x 0.5° at (7,6) w (-6.0° x -14.5°). For 
both fields, we fix the reddening at E(B — V) = 0.11, 
which is the mea n value over the core field, as judge d 
from the maps of iSchlegel. Finkbeiner. fc Davis! $1998). 
The left panel of Figure [5] shows the CMD of the Sgr 
core. This is a Hess difference that has been constructed 
by subtracting from the core Hess diagram the control 
field Hess diagram. To scale the control field density, we 
minimize the residual counts in the dashed box shown in 
Figure El There are a number of readily identifiable pop- 
ulations in the CMD, including the main sequence (MS), 
red giant branch (RGB), red clump (RC), red horizontal 
branch (RHB), blue horizontal branch (BHB), blue strag- 
glers (BS) and the luminous red giant branch (LRGB), 
together with M Giants selected in the same manner as 
in the right-hand panel of Figure [TJ We use the CMD 
boxes shown to build the LF of the core and calculate 
luminosity density profiles along the stream. It is reas- 
suri ng that our core CMD looks very similar to Figure 
1 of IBellazzini et ail (|2006l ). which has been constructed 
from the same data, but with a different algorithm. 

Figure [6] shows the observed LF of the Sgr remnant 
as derived from integrating the light in the CMD boxes. 
This drops dramatically beyond My ~ 5, and a number 
of possible extrapolations at the faint end are shown. 
However, mo st of the lum i nosity is in stars brighter than 
the turn-off. IMateo et~a l. (1995) computed the luminos- 
ity of the Sgr remnant by integrating the LF between the 
red clump and the upper-main sequence to get a surface 
brightness of 25.8 mag arcsec -2 . They account for the 
missing flux in fainter stars by extrapolating their LF 
to My — 13.2 to get 0.5 mag correction. They do not 
explicitly account for stars brighter than the red clump 
or bluer than the main-sequence turn-off. Repeating the 
calculation for our LF derived using the Bellazzini et al. 
(2006) data, we get a similar surface brightness in stars 
between the red clump and the upper main sequence, 
namely 25.7 mag arcsec -2 . However, we also find the 
BHBs, BSs, LRGBs and M giants contribute an addi- 
tional 0.3 mag. Two possible extrapolations at the faint 
end are shown, which contribute 0.5 mag (dot-dashed) 
and 0.15 (dotted) respectively. This can be compared to 
LFs of a typical globular cluster (red dashed, Cox 1999), 
the dSph Ursa Minor (black dashed, Wyse et al. 2002) 
and one of the most luminous globular clusters, M15 (red 
solid, Piotto et al. 1997). Taking our lead from the dSph 
LF suggests adding 0.5 mag for faint stars to arrive at a 
surface brightness of 24.9 mag arcsec -2 (compared to a 
value of 25.3 given in Mateo et al. 1995). 

3.3. The Density Profile of the Sagittarius Stream 

The CMD of streams A and B of t he Sgr lead i ng arm 
is already shown in Figure 4 of FBclokuro v et aTl ((2006a) 
and is known to look similar to that of the core. Here, we 
refine the calculation of the leading arm CMD in three 
ways. First, we use the late st data release DR7 of the 
SDSS (jAbazai lan eTaDHQOl). Second, we use a different 
comparison field, defined by reflecting the boundaries of 
the stream across the £ = 180° line (see Figure [2|). This 
assumes that the Galactic field star population is sym- 
metric across I = 180° and that the control fields there- 
fore probe a similar field star population as is present 
around the Sgr stream. Third, in order to combine por- 
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B-V B-V 

Fig. 5. — Comparison of the Hess diagrams of the Sgr core (left) and stream (right). Note that the populations in the core and the 
stream are very similar. Boxes are overlaid to identify characteristic populations used later in the paper, including main sequence (MS), 
red giant branch (RGB), red clump (RC), red horizontal branch (RHB), blue horizontal branch (BHB), blue stragglers (BS), and luminous 
red giant branch (LRGB). The black points in the left hand panel are M giant stars selected in the core using the cuts described above. 
The dashed box shown in the left panel is used to scale up the density in the control field. Ovcrplotted in the left panel are Dartmouth 
isochrones (Dottcr ct al. 2008) of varying age and metallicity (solid 15 Gyr, [Fe/H]= 0.0; dashed 11 Gyr, [Fe/H]=— 0.5; dash-dotted 8 Gyr, 
[Fe/H]= —1.0). We are confident that there are no turn-offs of older stellar populations in the Sagittarius core that were not detected in 
the Bellazzini et al. data. 



TABLE 1 

Heliocentric Distances to parts of the Sgr stream. SDSS and 2MASS observations are summarized in the top and bottom 

HALVES RESPECTIVELY. 



RA - Stream A 
Distance (kpc) 


215 
17 


210 
47 


205 
46 


200 
43 


195 
40 


190 
38 


185 
35 


180 
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30 
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27 
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24 
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RA - Stream B 
Distance (kpc) 
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RA - Trailing 
Distance (kpc) 
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tions of the stream at different distances, we offset each 
star in the stream and control field using a smooth map- 
ping between its right ascension and distance (see Figure 
[3]). Rather than showing the difference of two Hess dia- 
grams, here we choose to show the ratio. The advantage 
of this is that it accounts naturally for the fact that the 
depth of the CMD varies along the stream. This proce- 
dure helps to enhance fainter features like the horizontal 
branch and luminous part of the red giant branch, as can 
be seen in the right panel of Figure [5j This also explains 
the differences in relative density between the two pan- 



els. There are some small but significant differences in 
the location of CMD features. For example, the stream's 
CMD is slightly bluer than the core's. This could be due 
to metallicity gradients, although the bulk of it is prob- 
ably due to slight underestimation of extinction in the 
core (which is of course very difficult to measure in low 
latitude fields). Nonetheless, the size of our CMD boxes 
means that regardless of the slight color mismatch, Sgr 
member stars will still be selected. 

We count the numbers of stream stars in the popula- 
tion boxes marked in Figure [5j and subtract the counts 
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TABLE 2 

Single exponential fits to the A (top half) and B (bottom half) 

STREAMS (SEE FIGURE 0. THE LEFT COLUMN GIVES THE SLOPE OF THE 
EXPONENTIAL TOGETHER WITH ITS UNCERTAINTY. THE MIDDLE COLUMN 
GrVES THE LUMINOSITY PER KPC OF THE FITS AT RA= 220° . THE RIGHT 
COLUMN INDICATES THE NUMBER OF DATA POINTS USED IN THE 
EXPONENTIAL FIT. WHERE APPLICABLE THE VALUES IN BRACKETS 
INDICATE THE PARAMETERS DERIVED USING ALL DATA POINTS. 



Sample 




Exponent (xl0~ 3 ) 


— 

Luminosity (Lq x 10 ) 


Data Points 


Total 


7.0 ± 0.2 


27.00 


12 


MS 


6.6 ± 0.3 


15.92 


12 


BS 


4.5 ± 0.8 


0.24 


12 


RGB 


6.4 ± 0.5 


2.01 


12 


RC 


7.7 (3.5)± 1.0 (0.4) 


2.14 (1.66) 


6 


RHB 


8.1 ± 0.6 


3.18 


12 


BHB 


6.2 ± 1.0 


0.33 


12 


LRGB 


8.7 (-2.7)± 0.7 (0.5) 


2.12 (1.1) 


1 


Total 


4.5 ± 0.6 


15.01 


12 


MS 


4.9 ± 0.6 


9.20 


12 


BS 


2.7 ± 0.7 


0.18 


12 


RGB 


4.1 ± 0.7 


1.47 


12 


RC 


4.4 (0.5) ± 1.2 (0.4) 


1.2 (0.92) 





RHB 


5.1 ± 1.0 


1.47 


12 


BHB 


4.5 ± 1.2 


0.24 


12 


LRGB 


4.4 (-1.0) ± 1.2 (0.5) 


1.2 (1.2) 






each field's LF down to My = 13.2 by requiring that 
the slope of the stream's LF at My = 3.4 matches that 
of the core's extrapolated LF. This value is chosen as it 
corresponds to the limiting SDSS magnitude V ~ 22 in 
the most distant leading arm field. 

The Hess difference is integrated over each CMD box 
to give the total luminosity. Since each field has the 
same angular extent but is at a different heliocentric dis- 
tance, the luminosity is normalized by the length of the 
segment in kpc. The two panels of Figure [7] show the 
luminosity profiles as a function of right ascension for 
multiple populations in streams A and B, together with 
exponential fits. The choice of fitting law is motivated 
by the behavior of the total luminosity. We expect that 
all stellar populations behave in a similar way, and this is 
largely confirmed by our results. There are populations 
not so well described by an exponential law, but these 
deviations may be caused by noise in the data. Each 
of our fields probes at least some way down the main- 
sequence, so the only population that is affected by our 
LF correction is the main-sequence sample itself. How- 
ever, as shown in Figure [3 the RGB, RHB, BHB and 
BSs all have slopes that match the slope in the corrected 
main sequence sample, which suggests that the observed 
trend is not due to the correction we apply. Notice that 
the slopes of all the stellar populations in both streams 
A and B are positive at a statistically significant level. 
This implies that the corresponding densities of the pop- 
ulations are all fading with decreasing right ascension, 
that is increasing distance from the progenitor (as also 
noticed by Yanny et al. 2009). This is consistent with 
both A and B being part of the same leading stream. 

Stream A is roughly a factor two brighter in luminosity 
than B. About 60% (50%) of the signal is contributed by 
the main sequence stars when the dot-dashed (dotted) 
extrapolation of Figure [6] is used. The total luminosi- 
ties in the fields are on average a factor of 1.3 brighter 
when using the steeper extrapolation. BHBs and BSs 
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Fig. 6. — LF of the Sgr remnant (black), together with two pos- 
sible extrapolations at the faint end (dashed and dotted). Also 
shown are the LFs of Ursa Minor (black dashed), the globular 
cluster M15 (red solid), together with a mean globular cluster LF 
(red dashed). 



in the comparison field. Given the distance to each field, 
this produces a Hess difference in absolute magnitude 
versus color, My versus B — V. Due to the substantial 
distance gradient along the stream, the faintest popula- 
tion box - main sequence stars - is probed to different 
limiting magnitudes. To correct this, we make use of the 
fact that the stream's population is similar to that of the 
Sgr core, at least below the subgiant branch. We extend 
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Fig. 7. — Left and right panels show the luminosity per unit length in a variety of stellar populations in streams A and B respectively. 
The luminosities for the RGB, RC, RHB, BHB, BS and LRGB stars are derived directly from the Hess difference diagrams. The numbers 
are normalized to 1 kpc distance along the stream. The luminosities in the MS stars are corrected for missing flux using the derived LF. 
Error bars are determined assuming Poisson noise. A single exponential is fit to each population. However, the RC and LRGB populations 
show significant noise and hence for stream A only the first few data points follow an exponential fit. fn stream B the noise dominates 
these two populations and the exponential shown is estimated without using any data points as described in the main text. The fits to 
these two populations are therefore shown as dashed lines. The data points used in determining the fits are shown as filled symbols and 
those excluded as open symbols. The behavior of streams A and B is very similar, suggesting that both A and B may be parts of the same 
leading arm. Finally, we list at the side of the figure the mean of the fractional contribution to the total luminosity across the stream for 
each population. 



contribute little to the total luminosity, but can be very 
cleanly selected and so give a good indication of the un- 
derlying trend. By contrast, the counts corresponding 
to the RC and LRGB stars are much noisier. This is 
because the foreground contamination by the thick disc 
and halo is most severe, and varies considerably in the 
bright part of the CMD. In fact, the exponential fall-off 
is only present in the first few bins of stream A, corre- 
sponding to right ascensions between 180° and 205°. The 
dashed lines for stream A therefore show exponential fits 
using only the data at higher right ascension (shown as 
filled symbols) . For the fainter and therefore even noisier 
stream B, even this option is unavailable and so we as- 
sume a constant ratio of BHBs to RCs and LRGBs in the 
two streams to derive the dashed lines in the right-hand 
panel. This assumes that streams A and B have similar 
metallicities, which is supported by the work of Yanny et 
al. (2009). Table [5] gives the slope of each exponential fit 
together with its uncertainty, as well as the luminosity 
for each component at a right ascension of 220°. We also 
indicate the number of data points used in determining 
the exponential fit 

To obtain as complete a profile of the stream as pos- 
sible, we will need to exploit the complementarity of the 
2MASS and SDSS datasets. 2MASS pr ovides a power- 
ful pr obe of the trailing arm, for which Maicwsk i et al.l 
(2003) have already deduced a profile. There are also 
parts of the leading arm in 2MASS, for which we build 
the profile as illustrated in Figure [5J The upper panel 
shows the locations of M giants, selected using the same 



method as the zoom- in of Figure [JJ after offsetting stars 
to the distance of the core using Table [TJ The results 
are shown in the heliocentric coordinate system aligned 
with the Sgr orbital p lane, as defined by equation 9 of 
IMaiewski etaTl (|200l . There are two selection boxes 
in red to pick out the stream signal and the foreground. 
The resulting profiles of the raw number counts in 5° 
bins, the foreground and the corrected number counts 
are shown in the lower panel. Applying the same proce- 
dure to M giants in the trailing arm, we recover a profile 
similar to that of lMaiewski et akl (|2003|) , although rising 
in the outer parts. 

Before we can begin reassembly, we need to do three 
more things. First, we measure distance along the stream 
from the progenitor to the observed fields using Table [TJ 
Secondly, we combine streams A and B in the SDSS 
dataset into a single profile of the leading arm. Thirdly, 
in order to combine SDSS and 2MASS observations, we 
convert M giant number counts into total luminosity den- 
sities along the stream. We measure the relationship be- 
tween the number of M giants (-/V mo) and the tota l lumi - 
nosity at the locations where the IBellazzini et all (2006) 
and SDSS observations overlap with the 2MASS coverage 
(see the lower panel of Figure [3]). Here we can measure 
the total luminosity using the Hess difference method 
and count 2MASS M giants. The uncertainty in the ra- 
tio of the two is determined assuming Poisson statistics. 
There are clearly substantial gradients between the rem- 
nant and both the leading and trailing arms. This is 
illustrated in Figure [HJ We use a straight line fit, which 
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Fig. 8. — Top: Location of 2MASS M giants in the orbital plane defined by the Sgr orbit; see equation 9 of Majewski et al. (2003) 
for the definition of latitude and longitude (Bq,Aq). The red box encloses the Sgr leading (left) and trailing (right) streams, whilst the 
foreground box is chosen to have the same area. Bottom: Number counts (JVmg) of M giants in the stream in bins of 5° X 20° as a function 
of longitude in the orbital plane. Crosses (pluses) are counts in the central (foreground) boxes. The curves are estimates of the underlying 
smooth foreground. The red circles and triangles are the resulting profiles. The profile of the trailing arm derived by Majewski et al. (2003) 
is plotted in grey in the bottom right panel. 



has a gradient of 280±30Lq kpc -1 per M giant. Alterna- 
tively, assuming a step function constant in the remnant 
and constant in the stream does not affect the following 
analysis substantially. 

Figure 1101 shows the resulting luminosity profile as a 
function of distance along the stream for both the lead- 
ing (circles) and trailing (triangles) arms. Points de- 
rived from 2MASS data are represented by open sym- 
bols, whereas those from SDSS by filled symbols. The 
profiles of the leading and trailing arms are different by 
a factor of ~ 3 at small distances from the remnant, 
primarily as a consequence of the fact that the trailing 
arm is much more extended. Notice too the bump in 
the leading arm at ~ 40 kpc distance along the stream. 
Such featur es are common i n simulations - such as in 
Figure 1 of lLaw et al.1 §005) - as a consequence of the 
pile-up of material stripped at disk-crossings. There is 
also an evident break in the profile of the trailing arm 
at ~ 10 kpc distance a l ong t he stream, similar to that 
seen by lMaiewski et all (|2003l ). It is natural to associate 
this break with the transition between the core and the 
streams. Even though the profiles have a complex struc- 
ture, it is useful to give a simple fit. For the trailing arm, 
we fit a broken exponential and find indices —0.207 and 
0.008 either side of a break radius at 9.6 kpc. For the 
leading arm, the break radius lies behind the disk and 
we only measure the profile after the break as —0.007. 



3.4. The Luminosity of Sgr 

Now we begin the reassembly by calculating the to- 
tal luminosity in the Sgr remnant and tidal debris. We 
estimate the remnant's luminosity by summing the ob- 
servations of the trailing arm within the profile break to 
get 2.0 ±0.3 x 10 7 L Q . This is doubled on the assumption 
of the symmetric structure of the remnant, giving us a 
total magnitu de of My = — 14.2± 0.1, which is within 2 a 
of the value of iMateo et ail (fl998l) of My = -13.4 ±0.5. 

For the debris, we give two estimates. The first is a 
lower limit set by the data alone, namely 2.6±0.1 x 10 7 L Q 
for the leading arm and 8.2 ±0.7 x 10 6 L Q for the trailing. 
The combined value already contains more than half the 
core's luminosity, and is well constrained by the deep 
SDSS coverage of the leading arm. 

The second is an estimate of the luminosity contained 
within the entire debris stream, not simply the portions 
we observe. This is based on integrating the leading arm 
profile from 25 kpc (corresponding to the first 2MASS 
observation) to 300 kpc, which is triple the length to 
the last SDSS observation. This correspond s to a full 
wrap , as judged by the sim ulations (see e.g.. lLaw et al.l 
l200a IFellhauer et aLll200fih . This gives a luminosity of 
3.5±0.1xl0 7 L Q . This nominal error is so small because 
the fit is constrained by the SDSS observations which 
have small uncertainties due to the large number of stars 
in the SDSS sample. This also ensures that potentially 
larger photometric errors on individual stars are aver- 
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aged out. To account for the region between the core 
and 25 kpc, we average two possible extreme cases. If 
we assume that the core is symmetric, then the lead- 
ing arm profile needs to dip down sharply to meet the 
trailing arm profile at 10 kpc. This would contribute 
0.6 x 1O 7 L to the luminosity of the leading arm. Alter- 
natively, since we do not actually observe any indication 
of such a dip, the profile may continue rising until join- 
ing up with the core's profile at approximately 5 kpc. 
This would contribute 1.6 x 10 7 Lq. Therefore, for the 
missing inner portion of the leading arm, we take as a 
compromise the average value 1.1 ± 0.5 x 10 7 Lq where 
the large uncertainty reflects the range of possible values. 
The total leading arm luminosity (given by adding the 
extrapolated exponential and the estimate for the region 
between the core and the first observations)is then dou- 
bled to account for the trailing arm. We do not integrate 
the profile of the trailing arm since it is rising, which 
most likely reflects its approach to appocenter. Integrat- 
ing a rising profile is obviously very sensitive to the limits 
of integration. At present it is difficult to make a reason- 
able assumption about the behavior of the trailing arm 
far from the core. This finally gives for the total debris 
luminosity 9.2 ± 0.7 x 10 7 L©. This will be reduced to 
6.6 ± 0.7 x 10 7 L Q if the fainter LF in Figure[B]is used. It 
is evident that the formal error from Poisson noise, typi- 
cally of order of 10%, is much smaller than uncertainties 
induced by extrapolations. To account for the missing 
data, we were forced to extrapolate both the luminosity 
function and parts of the density profiles. 

The various luminosities used in the final answer are 
summarized in Table [3] The total luminosity of Sgr pro- 
genitor is 13.2 x lO 7 ^© or My ~ —15.5 with the brighter 
LF. If we repeat the above calculation with the fainter 
LF, we arrive at 9.6 x 10 7 L© or M v ~ -15.1. This im- 
plies that the progenitor is fainter than the SMC, but 
comparable to the brighter M31 dSphs like NGC 147 
(L = 9.4 ± 2.0 x 10 7 L Q ,M y = -15.1) and NGC 185 
(L = 13.5 ± 2.0 x 10 7 £ s , M v = -15.5). 

The systematic uncertainties clearly outweigh the Pois- 
son noise but we are confident that we have atleast found 
a firm lower limit (9.6 x 10 7 L©) to the luminosity. The 
sources of systematic uncertainty include the interpola- 
tion of the inner slope of the leading arm. This is over 
quite a modest range and unlikely to be seriously in er- 
ror. In addition considerable uncertainty derives from 
the extrapolation of the LF as shown above. Perhaps 
the greatest unknown is the behavior of the profiles at 
large distances. We have assumed that the exponential 
behavior continues to 300 kpc, but there may be pro- 
cesses that cause additional clumps and bumps. 

4. DISCUSSION AND CONCLUSIONS 

Using SDSS and 2MASS data, we have provided for 
the first time the luminosity profiles of both the leading 
and trailing arms of the Sgr stream. Prior to this work, 
the only available profile was that of the trailin g arm 
derived from M giants bv iMajewski et al.l (|2003f ). This 
was essentially flat, and so cannot be extrapolated to de- 
rive useful constraints on the mass lost during disruption. 
Here, our profiles for both the leading and trailing arms 
have gradients. It is the gradient in the leading arm that 
enables us to derive useful constraints. 

The luminosity in the leading stream clearly declines 
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Fig. 9. — Conversion factor between total luminosity and num- 
bers of M giants at different locations along the stream. The dash- 
dotted line shows the change in the conversion factor along the 
stream assuming that it is constant within the core, and constant 
along the stream. The dashed line shows the conversion factor 
assuming it varies linearly throughout the core and stream. 
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Fig. 10. — Comparison between the leading and trailing arms as 
seen in 2MASS and SDSS. We plot here the total luminosity per 
unit kpc along the stream versus the distance along the stream. 
For the 2MASS data, this is obtained from the number counts of 
M giants by the scaling described in the main text. The symbols 
have the same meaning as in the lower panel of Figure [3] For the 
trailing stream, the dashed line shows the best broken exponential 
fit to the data. For the leading stream, the slope of the inner part 
is fixed which determines the break point. 
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TABLE 3 

Luminosity of various parts of the Sgr stream together with 
their associated poisson uncertainties. 





Luminosity (in 


1O 7 L ) 


Stream Portion 


Faint LF 


Bright LF 


Core (Bellazzini et al.) 


2 x 1.5 ±0.3 


2 X 2.0 ±0.3 


Leading Arm (SDSS) 


2.0 ±0.1 


2.6 ±0.1 


Trailing Arm (2MASS) 


0.6 ±0.1 


0.8 ±0.1 


Leading Arm (25 to 300 kpc) 


2.5 ±0.1 


3.5 ±0.1 


Leading Arm (Core to 25 kpc) 


0.8 ±0.5 


1.1 ± 0.5 



with increasing distance from the progenitor. Although 
the beginning of the tail is absent because it lies behind 
the disk, the profile still shows mild evidence for a bump. 
This is most likely a pile-up at apogalacticon. The trail- 
ing stream can be traced all the way from the remnant 
to the anti-Center, and shows a clear break around the 
tidal radius as we move from the core to the tail. At first 
glance, this looks very different to the leading arm pro- 
file. But in reality it may just be a stretched-out version 
of the leading profile, as the trailing arm extends to much 
greater Galactocentric distances. In this picture, the ris- 
ing part of the profile is interpreted as the approach to 
an apogalacticon pile-up. We speculate that a similar 
but steeper rise would be observed in the leading arm, 
making the core profile symmetric, were it not obscured 
by the disk. 

It is worth emphasizing that a firm lower limit My ~ 
— 15 for the Sgr progenitor luminosity is given by the 
data with minimal assumptions. The most important 
is that the mass is equally apportioned between leading 
and trailing in both the remnant and the tails. Thanks 
to deep and wide SDSS coverage, the leading profile is 
measured with excellent accuracy. It is the tightness of 
the datapoints in the leading arm that is providing this 
firm lower limit. For example, integrating the exponen- 
tial profile out to 200 or 300 kpc makes little difference 
to the final answer. 

Our results may have implications for the scenario pro- 
posed by Fellhauer et al. (2006), who argued that stream 
A is material in the leading arm that was stripped re- 
cently, whereas stream B is debris in the trailing arm that 
was stripped long ago. Here, we have explicitly assumed 
that streams A and B are both parts of the same leading 
arm. T his is supported by recent findings of lYannv et al.l 
(2009) , who noted that the stars in streams A and B have 
very similar metallicities and velocities. Also, if stream B 
is old and trailing, then its density (about one half of A) 
is larger than the young trailing arm (about a third of A) , 
which seems counter-intuitive. So, although Fellhauer et 
al's (2006) interpretation has not been ruled out, there 
are possible causes for concern. 

The luminosity profiles offer a new, and largely unex- 
plored, constraint on the disruption process. The fact 
that the leading arm profile is not flat, but decreases, di- 
rectly limits the stellar mass lost in the disruption. The 
relative numbers of counts in the leading and trailing 
arms can be used to constrain the geometry of the or- 
bit. The stars move along the tail with a drift velocity 
deter mined by the proper ties of the progenitor and the 
orbit (De hnen et al.l [20041 ) . Therefore, the time of on- 
set of destruction of the stellar component is in principle 



recoverable. 

Our estimate of the progenitor luminosity is 9.6— 13.2 x 

10 7 L Q or My 15.1 - -15.5 with roughly 70% of the 

light in the tidal tails. The nearest look-alikes in the 
Local Group are probably the M31 satellites NGC 147 
and NGC 185. These two dSphs bracket properties of 
Sgr in terms of number of globular clusters, metallic- 
ity, gas content a nd velocity dispersion (see e.g. Tables 
12.4 and 12.6 of Ivan den B ergh 2001). Their mass-to- 
light ratios are comparatively modest at ~ 5 (Gch a et al.l 
2009). If we accept the comparison, then this suggests 
that the mass within the luminous radius of the Sgr pro- 
genitor was ~ 1O 9 M . This is right at the lower mass 
end of the range of models considered bv lJiang fc Binnevl 
(2000) who investigated the range of progenitor proper- 
ties which could reproduce the structure of the core rem- 
nant observed today. 

We can use the estimate that 70% of Sgr's luminos- 
ity now resides in the tidal debris to draw conclusions 
about the properties of the progenitor's dark matter 
halo. Peharrubia et al. (2008, ab) have modelled the 
disruption of dwarf galaxies with baryons following a 
King profile embedded in a dark matter Navarro-Frenk- 
White (NFW) halo. The simulations suggest that the 
evolution of the parameters of the King profile is gov- 
erned primarily by the amount of mass that has been 
stripped, rather than the details of the dwarf's orbit or 
the segregation of baryons and dark matter. Using as 
the present parameters of the bound core those deter- 
mined by Majewski et al. (2003), we find that prior to 
disruption the King profile had a central velocity dis- 
persion (To = 23 kms -1 , projected central stellar density 
So = 5.7 x 1O 6 M kpc -2 and a core radius Rc = 1.5 
kpc. This implies that prior to mass stripping the cen- 
tral velocity dispersion, projected central stellar density 
and core radius were, respectively, a factor 2.0, 5.3, and 
1.2 higher than at present. Subsequently, we apply the 
method of Pefiarrubia et al. (2008a), to estimate the 
properties of the NFW dark matter halo in which Sgr 
was embedded prior to mass stripping implied by this 
King profile. These authors solve the Jeans equations 
in order to find a set of r Ci , nax — u c ,max parameters (the 
maximum circular velocity of the halo and the radius at 
which this is achieved which fully specify the NFW halo) 
that would allow a given velocity dispersion of the stellar 
component. This analyses yields a family of NFW mod- 
els depending on the assumed spatial segregation of the 
stellar component within the dark matter halo. To break 
this degeneracy the authors appeal to the results of cos- 
mological N-body simulations which show a strong corre- 
lation between an NFW halo's u c ,max and r cmax . Using 
the values for the King profile described above we de- 
rive a peak circular velocity U c ,max = 35 kms -1 and scale 
radius rs = 3.5 kpc ( 

^c,max ~ 2.17rg). The maximum 
circular velocity and scale radius are consistent with the 
values for the NFW haloes of other classical dSphs (see 
Figure 5 of Peharrubia et al. 2008a), but already towards 
the high mass end, with M vir « 1.0 x 1O 1O M . The Draco 
dSph, found by Peharrubia et al. to have the most mas- 
sive dark matter halo, has M v i r ps 6.3 x 1O 9 M . As 
shown by Peharrubia et al. (2006), at this mass, dynam- 
ical friction can reduce the orbits apo- and pericenters 
by a factor of > 2 over a Hubble time. We note that 
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Fig. 11. — Virial mass of the NFW halo at z = as a function 
of luminosity remaining in the core. The shaded region shows the 
range of luminosities determined in this paper. If we are missing 
light from additional wraps or apocenter pile-ups the mass quickly 
exceeds that of the other classical dwarfs. The dashed line indicates 
the virial mass of Draco (the most massive classical dwarf according 
to Penarrubia et al. 2008a). 

9.6 — 13.2 x 10 7 Lo is likely a lower bound for the progen- 
itor's luminosity since we have only considered one wrap 



of debris and have neglected the possibility of additional 
pile-ups at previous apocenters in the orbit. The amount 
of light contained in these will depend on the details of 
the mass loss process. With additional light in the de- 
bris Sagittarius quickly becomes more massive than any 
of the other known Milky Way dwarfs (see Figure ITT|) . 
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